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The GLAST sateUite mission wiU study the gamma ray sky with considerably greater exposure 
than its predecessor EGRET. In addition, it will be capable of measuring the arrival directions 
of gamma rays with much greater precision. These features each significantly enhance GLAST's 
potential for identifying gamma rays produced in the annihilations of dark matter particles. The 
combined use of spectral and angular information, however, is essential if the full sensitivity of 
GLAST to dark matter is to be exploited. In this paper, we discuss the separation of dark matter 
annihilation products from astrophysical backgrounds, focusing on the Galactic Center region, and 
perform a forecast for such an analysis. We consider both point-like and diffuse astrophysical back- 
grounds and model them using a point-spread-function for GLAST. While the results of our study 
depend on the specific characteristics of the dark matter signal and astrophysical backgrounds, we 
find that in many scenarios it is possible to successfully identify dark matter annihilation radiation, 
even in the presence of significant astrophysical backgrounds. 

PACS numbers: 95.35.-|-d; 95.85.Pw 
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I. INTRODUCTION 

Dark matter constitutes an essential element of mod- 
ern cosmology. Evidence of its existence come from 
a wide variety of observations including the rotational 
speeds of galaxies J]| , the orbital velocities of galax- 
ies within clusters [2], gravitational lensing Q, the cos- 
mic microwave background [j], the light element abun- 
dances [5] and large scale structure Q ■ But despite these 
many observational indications of dark matter, we re- 
main ignorant of the particle nature of this substance. 
To reveal the identity of dark matter, it will be cru- 
cial to measure its non-gravitational couplings. Efforts 
in this direction include direct detection experiments, 
which are designed to observe the elastic scattering of 
dark matter particles with nuclei, and indirect detec- 
tion experiments, which aim to detect the annihilation 
products of dark matter. Such annihilations could poten- 
tially produce observable fluxes of gamma rays, neutri- 
nos [3], positrons Q, antiprotons 9], antidcuterons \W^ . 
synchrotron radiation [ll[ and X-rays [l^ . Of these vari- 
ous channels, gamma rays have the important advantage 
of retaining directional information. This is a feature 
that can be used, together with the peculiar energy spec- 
tra expected from dark matter annihilations, to disentan- 
gle dark matter annihilation products from astrophysical 
backgrounds p^ . 

In this paper, we study the ability of gamma ray 
experiments to identify dark matter annihilation ra- 
diation from the Galactic Center region (for earlier 
work on this subject, see Refs. [Ijl) by using both 
spectral and angular information. This is a goal of 
existing Atmospheric Cerenkov Telescopes (ACTs) in- 
cluding HESS m, MAGIC tia], VERITAS [13 and 
CANGAROO-inlill as well as of the forthcoming satel- 
lite mission, GLAST [l^. Motivated by its sensitivity to 



gamma rays over an energy range well suited to dark mat- 
ter searches (0.1-300 GeV), we focus here on the GLAST 
experiment. 

The remainder of this paper is structured as follows: 
In Sees. |TT] and IIIII we describe the model used for the 
point-like and diffuse astrophysical backgrounds and the 
expected characteristics of a dark matter annihilation 
signal, respectively. In Sec. IIVI we describe our analy- 
sis method and assess the ability of GLAST to either 
place limits on or identify dark matter annihilations in 
the Galactic Center region. In Sec. |V] we extend our 
method to assess GLAST's ability to measure the an- 
nihilation cross section, mass and distribution of dark 
matter. In Sec. IVII we summarize and present our con- 
clusions. 



II. MODELING THE GALACTIC CENTER 
BACKGROUNDS 

The Galactic Center is a complex region of the sky at 
all wavelengths, the gamma-ray window being no excep- 
tion. In this section, we discuss how, in our analysis, 
we treat the backgrounds for dark matter searches due 
to known and unknown astrophysical sources of gamma 
rays. 

The first of the backgrounds we consider is the rel- 
atively bright, very high-energ y g amma ray source ob- 
served by HESS M, MAGIC [lij, WHIPPLE [H and 
CANGAROO-II [23]. This source is consistent with 
point-like emission and is located at / — 359°56'41.1" ± 
6.4" (stat.), b = -0°2'39.2"±5.9" (stat.) with a system- 
atic pointing error of 28" [2^ . It appears to be coincident 
with the position of Sgr A*, the black hole constituting 
the dynamical center of the Milky Way. The spectrum of 
this source is well described by a power-law with a spec- 



tral index of a = 2.25 ± 0.04(stat) ± O.lO(syst) over the 
range of approximately 160 GeV to 20 TeV. Although 
speculations were initially made that this source could 
be the product of annihilations of very heavy ('^10 to 
50 TeV) dark matter particles J25(], this interpretation is 
disfavored by the power-law form of the observed spec- 
trum. The source of these gamma rays is, instead, likely 
an astrophysical accelerator associated with our Galaxy's 
central supermassive black hole |26| . In our analysis, we 
treat this source as a background for dark matter searches 
(see also Ref. ^). 

Following the measurements of HESS, we describe the 
spectrum of this source as a power-law given by: 



-2.25 



$ACT^l_OxlO-M^J GeV-icm-2s-i. (1) 



At energies below ~ 200 GeV, the spectrum of this source 
has not yet been measured. GLAST, however, will be ca- 
pable of measuring the spectrum of this source at energies 
below the thresholds of HESS and other ACTs. As any 
signal from dark matter annihilations is expected to be 
spatially extended, as a result of the dark matter halo 
profile (see Sec. IHip . GLAST's spectral measurement in 
the angular window around the source could be used to 
obtain a relatively pure determination of the background 
spectrum from the HESS source. This approach would 
be less effective, however, in the case in which the spatial 
distribution of dark matter annihilations is highly con- 
centrated, such as might occur for a highly adiabatically 
contracted halo profile |23] or for a density spike resulting 
from adiabatic accretion of dark matter onto our galaxy's 
supermassive black hole [29|. In any case, the spectral in- 
dex of this source has been measured by HESS at energies 
well beyond those expected to be relevant to data matter 
searches. Unless the spectral index changes significantly 
over the energy range of interest to GLAST, this infor- 
mation can be used to effectively predict the spectrum 
GLAST will detect from this astrophysical source. 

In addition to the HESS source, a yet unidentified 
source has been detected by EGRET approximately 0.2° 
away from the dynamical center of our galaxy [30|, '3l|. 
Although the spectrum of this source is not yet well mea- 
sured, its spectral index appears to be similar to that of 
the HESS source {a « 2.2). We model the flux from this 
source as: 



$^^ = 2.2x10-^ i'Fv) GeV-^cm-^s-i 
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X exp(-£;^/30GeV). 



(2) 



EGRET detected the presence of this source at ener- 
gies up to ~ 10 GeV. If its power-law spectrum continues 
up to energies detectable by HESS and other ACTs, they 
would have observed the source as well. As they do not, 
we are forced to introduce a cutoff in the spectrum of 
this source. The precise value of the cutoff (30 GeV) was 
chosen arbitrarily. Once again, GLAST will measure the 
spectrum of this and other point sources and will thus 



not have to rely on such extrapolations or speculations 
in an actual analysis. 

Even ideal point-like sources appear to be slightly ex- 
tended in a realistic experiment due to the finite angular 
resolution of the detector. To account for this, we model 
the GLAST point spread function by: 
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where is the difference between the measured and actual 
directions of the observed gamma ray and the solid angle 
element is dfJ = 6d9d(j), 4> being the azimuthal direction 
of the displacement with respect to an arbitrary azimuth. 
The function 5{E^) is the angle within which 68% of the 
gamma rays are reconstructed and is given in degrees by: 

\og^^[5{E^)] « -0.6 - 0.81ogio[i?^/l GeV]. (4) 

Strictly speaking, the normalization factor in Eq. ([3]) is 
correct only in the planar limit sinS ^ 6, that is when 
(5 ^ TT, which for the case at hand is an excellent approx- 
imation, see Eq. (H]). In this limit, the integral of Eq. ([3]) 
over the whole solid angle is equal to unity. 

Finally, we consider the diffuse background of gamma 
rays in the region surrounding the Galactic Center. Un- 
like in the case of point sources, the detailed angular dis- 
tribution of the diffuse background is not known and thus 
cannot be easily separated from a potential dark matter 
annihilation signal using angular information alone. In 
our analysis we consider two extreme possibilities in this 
regard. In the first case, we use a flat distribution for 
the diffuse background, with the same intensity in each 
angular bin (over a 2° x 2° window centered around the 
Galactic Center). In the second case, we use an angular 
distribution which matches that of the dark matter sig- 
nal. The latter case represents the extreme limit where 
there is no angular information with which to separate 
the annihilation signal from the diffuse background. Al- 
though the most realistic situation is bracketed by these 
two limits, we expect the uniform background case to be 
closer to reality than the latter one. 

The spectrum of the diffuse background is also not 
known in advance, and may be difficult to separate from 
dark matter annihilation radiation. Observations by 
HESS of the diffuse emission from the Galactic Center 
Ridge [321 find a spectrum in the region of the sky given 
by -0.8° < I < 0.8°, -0.3° < 5 < 0.3° which can be 
described by a power-law: 
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The HESS observations do not, however, constrain the 
diffuse spectrum (or angular distribution) at energies be- 
low '^200 GeV, which arc of the most interest to GLAST. 
In our analysis, we assume the general parameterization 
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where a is allowed to vary between 1.5 and 3.0. We adopt 
an overall normalization, A, such that the integrated flux 
of the diffuse background between 1 GeV and 300 GeV 
in a 2° X 2° field of view around the Galactic Center is 
equal to 10"'* cni~^ s~* sr~*. Wc do not, however, assume 
that this normalization is known in our analysis, leaving 
open the possibility that some of the diffuse gamma rays 
observed are the product of dark matter annihilations. 



zenith gives a mean effective area equal to 60% of the 
peak area [33|. We quantify the latter two effects by 
multiplying A^s by the exposure parameter e, with e = 
0.5 X 0.6 ~ 0.3. The two distinctive features seen in the 
figure correspond to the HESS and EGRET sources de- 
scribed above. Also included is a diffuse background, dis- 
tributed isotropically, with a spectrum given by $'^'^ = 
1.4 X 10-'^{E^/GeV)-^-^ GeV"* cm'^ s"* sr"*. 
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FIG. 1: A simulated sky map of the gamma ray backgrounds 
present in the region of the Galactic Center, after two years of 
observation by GLAST. Each point denotes one gamma ray 
detected. In the top and bottom frames, photons with energy 
above 3 and 10 GeV are shown, respectively. 

In Fig. [1] we show an example of the simulated back- 
grounds (including Poisson noise) as might be observed 
by GLAST in two years of observation. For GLAST, 
we have used a constant effective peak area above 1 
GeV of Acs = 8.5 x 10'^ cm^. We have also assumed 
that the Galactic Center will be within GLAST ficld-of- 
view 50% of the time and that the reduction in effective 
area for sources which are not located at the instrument 



III. GAMMA RAYS FROM DARK MATTER 
ANNIHILATIONS 

The energy and angular dependent flux of gamma rays 
produced in dark matter annihilations is described by 



$ 



DM 



diV^ {av) 



p^{r)dl, 



(7) 
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where {av) is the WIMP annihilation cross section multi- 
plied by the relative velocity of the two WIMPs (averaged 
over the WIMP velocity distribution) , mx is the mass of 
the WIMP, -ip is the angle observed relative to the direc- 
tion of the Galactic Center, p{r) is the dark matter den- 
sity as a function of distance to the Galactic Center, and 
the integral is performed over the line-of-sight. dNj/dEj 
is the gamma ray spectrum generated per WIMP annihi- 
lation. The spectrum of gamma rays produced through 
dark matter annihilations depends on the nature of the 
WIMP. Neutralinos, for example, typically annihilate to 
final states consisting of heavy fermions (65, it, t~^t~) 
or gauge and/or Higgs bosons {ZZ, W^W^ , HA, hA, 
ZH, Zh, ZA, W^H^, where H, h, A and iJ± are the 
Higgs bosons of the Minimal Supersymmetric Standard 
Model) [3J|- With the exception of the t^t^ channel, 
each of these annihilation modes result in a very similar 
spectrum of gamma rays. In Fig. [21 we show the pre- 
dicted gamma ray spectrum, per annihilation, for several 
possible WIMP annihilation modes. In this article, we 
do not consider mono-energetic gamma ray lines [35|, as 
they are expected to produce far fewer events in GLAST 
than continuum emission. 

In convenient units, Eq. ([7]) can be recast as: 



$^(£:^, V) ~ 2.8 X 10-*'^cm-^s 



dN^ 



(av) 



100 GeV 



X J(V') (8) 



dE^ \3 X IQ-^^cm^/s J \ mx 

where the dimensionless function J{ip) depends only on 
the dark matter distribution in the halo and is defined 
by convention as: 
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(9) 

To calculate J(V'), a model for the dark matter halo dis- 
tribution must be adopted. A commonly used parame- 
terization of halo profiles is given by 



p(r) = 



(r/i?)''[l + (r/i?)"](/5-7)/c 



(10) 



t3 



loi t^ 



10^ 



^ 10-1 



10" 



I ' ' ' ' I r 




mx=100 GeV 

I , , I .... I L. 



5 10 20 50 100 

E^ (GeV) 



-■:--'.:._; . ' I ' ' "I 

10^ r - ^ '^^ 'i»,. 



I ' ' "I 



^,, bb, tt, V*V', ZZ, HA, ZH 




FIG. 2: The gamma ray spectrum per annihilation for a 100 
GeV (top) and 500 GeV (bottom) WIMP. Each curve denotes 
the result for a different dominant annihilation mode. 



where i? ~ 20 kpc is the scale radius and po is fixed by im- 
posing that the dark matter density at the distance of the 
Sun from the Galactic Center is equal to 0.3 GeV/cm'^. 
Among the most frequently used models is the Navarro- 
Frcnk- White (NFW) profile, which is described by a = 1, 
/3 = 3 and 7=1 [36]. When considering the region of 
the Galactic Center, the most important feature of the 
halo profile is the inner slope, 7, which can be steeper or 
harder than the assumed fiducial value 7 = 1. For ex- 
ample, the Moore et at profile is described by a = 1.5, 
/3 = 3, 7 = 1.5 [13. Note that for 7 > 1.5, the integral 
in Eq. ([9]) diverges. To avoid this behavior in the case 
of the Moore et al. profile, we impose a flat core within 
10~^ pc of the Galactic Center. 

Although we will use the NFW and Moore et al. halo 
profiles as benchmarks, they certainly do not exhaust all 



possibilities. For a number of reasons, it is quite diffi- 
cult to predict the dark matter distribution in the inner 
parsecs of the Galaxy, in which we are most interested. 
The resolution of N-body simulations is limited to scales 
of approximately ^ 10^ parsecs or so. Furthermore, the 
gravitational potential in the inner region of the Milky 
Way is dominated not by dark matter, but by baryons, 
whose effects are not included in such simulations. The 
precise impact of baryons on the dark matter distribu- 
tion is difficult to predict, although an enhancement in 
the dark matter annihilation rate due to adiabatic com- 
pression is generally expected [28|]. The adiabatic accre- 
tion of dark matter onto the central supermassive black 
hole may also lead to the formation of a density spike in 
the dark matter distribution. Such a spike would likely 
result in a very high dark matter annihilation rate '29!| . 
The remaining factors in Eq. ([8]) depend on the parti- 
cle physics of the dark matter candidate, i.e. the mass, 
cross section, and spectrum per annihilation. As a bench- 
mark value, we adopt a WIMP annihilation cross section 
of {av) « 3 X 10~^^ cm'^/s. This is motivated by the 
fact that a WIMP annihilating with such a cross section 
during the freeze-out epoch will be generated as a ther- 
mal relic with a density similar to the measured dark 
matter abundance [^. If these annihilations take place 
largely through S'-wave processes, then the annihilation 
cross section of WIMPs in the Galactic halo (ie. in the 
low velocity limit) will also be approximately equal to 
this value. The annihilation cross section could be con- 
siderably reduced, however, if P-wave processes or coan- 
nihilations between the WIMP and other particles are 
significant during the freeze-out process. 



IV. PROJECTED CONSTRAINTS 

In order to evaluate the ability of GLAST to identify 
or constrain the properties of annihilating dark matter, 
we consider the spatial and spectral information simul- 
taneously. We use the model introduced in the previous 
sections to construct a simulated sky map, free of a dark 
matter signal. Then, for given dark matter mass and 
annihilation cross section, we compute the x^ of the sim- 
ulated sky against models including a contribution from 
dark matter annihilation radiation. We then repeat the 
analysis against the same simulated sky for a wide range 
of models. 

The simulated sky contains the two resolved sources 
and the diffuse background described in Section [ill For 
a given angular pixel, v, and energy bin, k, we calculate 
the total number of expected photons: 
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(11) 



and then add in Poisson noise. For a given running time 
of the experiment, i, each of these terms is obtained by 
integrating the flux over an given energy and angular bin, 

dV 



B%=fKex 



dn 



An„ 



dSAeff$«^. (12) 



Each angular bin covers solid angle A^^, = 0.01° x 0.01° 
(smaller than the angular resolution of GLAST). Since we 
consider the inner 2° x 2° field-of-view, there are a total 
of 40,000 bins. We include 25 energy bins distributed log- 
arithmically between 1 and 300 GeV (approximately 10 
bins per energy decade, which is realistic given GLAST 
energy resolution). 

The distribution over the sky, 0°^^ , therefore, does not 
contain any contribution from dark matter annihilations. 
We compare this to a model which includes both the 
backgrounds and a signal from dark matter: 

R.,k = 5.,fe(mx, {av)) + B^2^ + B^^ + ^'fe (^'- «') • 

(13) 
We do this by computing the x^ of the fit as follows ^ : 
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(14) 

The free parameters in the y^ are the slope of the dif- 
fuse background, the amplitudes of the point source back- 
grounds, and the mass and annihilation cross section of 
the dark matter particle. We fix the amplitude of the 
diffuse spectrum by summing over all events. A more 
complete analysis would marginalize over this amplitude 
as well, but we are in the limit where the overall number 
of photons is dominated by the astrophysical background. 
Fixing the quantity reduces the CPU time required while 
leading to only a small underestimate of the errors on the 
dark matter parameters. 

Although we have included only the HESS and 
EGRET sources (in addition to the diffuse background) 
in our analysis, this method could easily be extended to 
include any other astrophysical point sources to be dis- 
covered by GLAST. Also, in principle one could have 
introduced additional parameters in the modeling of the 
point-like sources. Computational constraints, however, 
require us to limit the number of free parameters in our 
analysis, and the most difficult background to separate is 
the diffuse component. Concerning the point-like source 
at the Galactic Center, it is worth mentioning that the 
bulk of the statistical significance of the dark matter 
annihilation signal does not come from the inner 0.1° 
around the Galactic Center (where the HESS source dom- 
inates), but rather from the surrounding angular region, 
even in the case of a cusped halo profile ^] . 

By scanning over this five-dimensional parameter 
space, we can evaluate how large a contribution from dark 



^ Strictly speaking, this expression is valid only if several events 
are present in a given energy/angular bin. This is not always 
the case in the present application. The most important bins, 
however, which contribute most to the discriniating power of this 
calculation, do contain a sufficient number of events. We have 
checked this by increasing the size of our angular bins and find 
good agreement. 



10" 



-24 



m 


10-26 






> 

b 


10-28 



10" 



-30 



1 
NFW 




1 , , , , 


1 1 ' ' 


1 


zsss^ 


- 








- 




. ' "-li 




- 






1 !. , , , 




1 , , , , 



10 20 50 100 200 500 1000 
mx (GeV) 





10-24 


- 


M 


1 
oore 


' 1 ' ' ' ' 1 1 i i 1 i i i i 
et al. 


Ui 


10-26 


_ 








'^^^^^^yi. ■:,■ :-. 


> 

b 


10-28 

in-30 


- 


- 


— 


-— 


, 1 !. , , , 1 1 , , 1 , , , , 



10 20 



50 100 200 
mx (GeV) 



500 1000 



FIG. 3: The projected exclusion limits at 95% confidence level 
from GLAST (after ten years) on the WIMP annihilation 
cross section, as a function of the WIMP mass. The region 
above the dotted line is already excluded by EGRET 30]. 
The dashed and solid lines show the projections for GLAST 
for an assumed isotropic diffuse background and the limit case 
where the astrophysical background has exactly the same an- 
gular distribution of the DM signal, respectively. In the upper 
and lower frames, the NFW and Moore et al. halo profiles have 
been adopted, respectively. Also shown are points represent- 
ing a random scan of supersymmetric models. 



matter annihilations can be contained in Ri,^k before the 
total signal becomes inconsistent with the background- 
only sky, 0°*^ . In Fig. [H we show the projected exclu- 
sion limits at the 95% confidence level in the {mx, {cv)} 
plane for ten years of collection time by GLAST and for 
the cases of an NFW (top frame) or Moore et al. (bottom 
frame) halo profile and WIMP annihilations to bb. The 
solid line in each frame represents the limit found if the 
diffuse background is assumed to be distributed isotrop- 



ically, while the dashed hne represents the conservative 
hmit obtained if the diffuse background has the same an- 
gular distribution as the dark matter signal {i.e. the case 
in which angular information is not useful in disentan- 
gling the signal from the diffuse background). For values 
of (av) below the corresponding lines, a pure background 
model of the kind of Eq. (fTT|) is expected to be consistent 
with the data. The fact that the limits are significantly 
stronger in the uniform background case is the manifes- 
tation of the improved sensitivity which can be achieved 
by an analysis including both energy and angular infor- 
mation. For comparison, in Fig. [51 we also show the re- 
gion already excluded by EGRET [30| (above the dotted 
line) and the mass and cross section of neutralino models 
found in a random scan over supersymmetric parameters, 
as calculated using DarkSUSY 39]. As expected, many 
of the models cluster around {<tv) ^ 3 x 10^^^ cm'^/s, 
the value required of a thermal relic annihilating via an 
5-wave amplitude. Each point shown represents a model 
which respects all direct collider constraint and generates 
a thermal dark matter abundance consistent with the ob- 
served dark matter density. In our scan, we have varied 
the SUSY parameters Af2, |mI ^-nd rriq up to 2 TeV, uia 
and nil up to 1 TeV and tan (3 up to 60. We have assumed 
the gaugino masses evolve to a single unified scale, such 
that Ml « 0.5M2, Mg w 2.7 M2. 



V. MEASURING THE PROPERTIES OF DARK 
MATTER 

Once gamma rays are identified as having been pro- 
duced in dark matter annihilations, such observations 
could then be used to measure the characteristics of 
the dark matter particle, including its mass, annihila- 
tion cross section and spatial distribution. Such deter- 
minations are an important step toward identifying the 
particle nature of dark matter. In this section, we discuss 
GLAST's ability to constrain these properties. 

To accomplish this, we do a similar calculation to that 
performed in Sec. lIV| but now also include a contribution 
from dark matter annihilations in the quantity, 0°^^. In 
particular, we include the contribution from WIMPs dis- 
tributed with an NFW halo profile, with an annihilation 
cross section of 3 x 10~^^ cm^/s and a mass of 100 GeV. 
We then calculate the statistical significance at which 
these properties can be measured by GLAST. 

In Fig. m we show the ability of GLAST to determine 
the WIMP mass and annihilation cross section in this 
case. In each frame, we show the projected 2 and 3 a 
constraints on these properties as may be determined by 
GLAST, assuming an isotropic diffuse background (in ad- 
dition to background point sources). In the top frame, we 
treat the shape of the halo profile (NFW) as if it is known 
in advance. Of course, this is not a realistic assump- 
tion, and a less accurate determination of the WIMP 
mass must be expected in a more realistic treatment. In 
the lower panel of Fig. 2] we report the results obtained 
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FIG. 4; The ability of GLAST to measure the annihilation 
cross section and mass of dark matter after ten years of ob- 
servation. Here, we have used a benchmark scenario with 
mx = 100 GeV, (av) = 3 x IQ-^'' cm^s and an NFW halo 
profile. The inner and outer contours in each frame represent 
the 2 and 3 a regions, respectively. In the top frame, the halo 
profile shape was treated as if it is known in advance. In the 
lower frame, we marginalize over the inner slope of the profile. 



marginalizing over the inner slope of the halo profile, 7 
[see Eq. p^ ]. Without a prior on the inner halo slope, 
the constraint on the dark matter mass worsens by a fac- 
tor -2. 

If the spectrum and angular distribution of gamma 
rays from dark matter annihilations in the Galactic Cen- 
ter region are sufficiently well measured, it will also be 
possible to measure the underlying dark matter distri- 
bution. In Fig. [5l we project the results from the lower 
frame of Fig. 3] in the {mx, 7} plane (marginalized over 
the annihilation cross section) . In our benchmark model, 
the inner slope of the halo profile can be determined at 
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FIG. 5: The ability of GLAST to measure the inner slope 
of the halo profile and the mass of dark matter particle 
(marginalizing over the annihilation cross section) after ten 
years of observation. Here, we have used a benchmark sce- 
nario with rax = 100 GeV, {pv) = 3 x 10^^^ cm'^/s and an 
NFW halo profile. The inner and outer contours represent 
the 2 and 3 a regions, respectively. 



approximately the ~ 10% level. 



VI. DISCUSSION AND CONCLUSIONS 

The most challenging task for indirect dark matter 
searches is not detecting the products of dark matter an- 
nihilations, but in confidently identifying those particles 
as such. In particular, any signal must be separated from 
astrophysical backgrounds if it is to be reliably claimed 
to be a detection of dark matter annihilation products. 
This is certainly true in the case of gamma ray telescopes 
hoping to observe dark matter annihilations in the region 
of the Galactic Center, where astrophysical backgrounds 
are especially foreboding. 

In this paper, we have studied quantitatively the abil- 
ity of GLAST to identify gamma rays from dark matter 
annihilations in the Galactic Center region, using all in- 
formation encoded in the spectrum and angular distribu- 
tion of the observed events. Due to the complex nature 
of the Galactic Center's backgrounds and dark matter 
distribution, it is important for GLAST to make full use 
of this information to best separate dark matter annihila- 
tion products from any astrophysical backgrounds which 
are present. 

When GLAST begins its mission in 2008, it will de- 
tect a number of astrophysical sources in the region of 



the sky around the Galactic Center. These include the 
point sources identified by HESS and EGRET, and per- 
haps other, thus far unknown, sources. A diffuse gamma 
ray background will also likely be present. In this arti- 
cle, we have modeled these backgrounds based on known 
spectral information, and using a point-spread-function 
for GLAST, in the hope of identifying statistically sig- 
nificant departures from this model resulting from dark 
matter annihilations. 

The spectrum of gamma rays from dark matter an- 
nihilations is quite distinctive, and does not resemble 
the power-law form observed from typically astrophysical 
sources. Furthermore, the angular distribution of dark 
matter annihilation products is not a simple point source, 
nor is it isotropic or trace the Galactic Ridge, as the back- 
grounds are expected to. By exploiting these differences, 
it may be possible to identify the products of dark mat- 
ter annihilation observed from the inner galaxy, even they 
are a subdominant component of the total emission. 

If dark matter annihilation products are identified by 
GLAST, then it may also be possible to measure or con- 
strain the properties of dark matter, including its mass, 
annihilation cross section and spatial distribution. We 
find that GLAST is unlikely to determine the WIMP's 
mass with high precision, however. For example, for the 
case of a 100 GeV WIMP with an annihilation cross sec- 
tion of 3 xlO~^^ cm'^/s and distributed with an NFW halo 
profile, the mass could be determined to lie within ap- 
proximately 50-300 GeV. In the same benchmark model, 
the inner slope of the dark matter halo profile could be 
determined to ^ 10% precision by GLAST. 

Note added: After we had completed this study, we 
became aware of updated estimates for the performance 
of GLAST 40]. In particular, the point spread function 
that we have used is a factor 3 -;- 4 more optimistic than 
was found in the current estimate. As a consequence, a 
reduction of the ability of GLAST to identify dark matter 
and to measure its properties is expected, but not to an 
overwhelming extent: we checked for example that for a 
NFW profile and a DM mass of 100 GeV the bound wors- 
ens by less than 30%. Also, additional effects not taken 
into account here (like effective area degradation due to 
orbit inclination and spacecraft rocking; non-negligible 
dead-time and a high-particle-background associated to 
orbital passage over the South-Atlantic Anomaly) may 
reduce the exposure by up to a factor ~ 2 compared to 
the present estimate. We would like to thank the referee 
for bringing this to our attention. 
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